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AERODYNAMIC INVESTIGATION OF A PARABOLIC BODY OF 


REVOLUTION AT MACH NUMBER OF 1.92 AND SOME EFFECTS OF 


AN ANNULAR JET EXHAUSTING FROM THE BASE 


By Eugene S. Love 


SUMMARY 


An aerodynamic investigation of a parabolic body of revolution 
was conducted at a Mach number of 1.92 with and without the effects 
of an annular jet exhausting from the base. Measurements with the jet 
off were made of lift, drag, pitching moment * radial and axial 
pressure distributions, and base pressures. With the jet in operation, 
measurements were made of the pressures over the rear of the body with 
the primary variables being angle of attack, ratio of jet velocity 
to free— stream velocity, and ratio of jet pressure to stream pressure*. 

The results with the Jet inoperative showed that the radial 
pressures over the body varied appreciably from the distribution gen- 
erally employed in approximate theories. The linearized solutions for 
lift, pitching moment, and center of pressure gave relatively poor 
predictions of the experimental results. An analysis of several 
theoretical methods for calculating pressure distribution and wave 
drag showed some methods to give results considerably in disagreement 
with experimental values. 

Maximum effects of the jet were obtained at the lower ratio of 
jet velocity to stream velocity and the highest ratio of jet pressure 
to stream pressure. These effects amounted to a slight decrease in 
fore drag, a reduction in lift, and a shift of center of pressure in 
a destabilizing direction. 
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INTRODUCTION 


Aerodynamic investigations at supersonic speeds of "bodies of 
revolution simulating those containing jet— propulsion units have almost 
entirely neglected the effects of the jet flow upon the flow over the 
rear of the "body. An experimental subsonic investigation of the effects 
of the jet upon the aerodynamic characteristics of the aggregate 
A— 5 missile ("body plus four equally spaced tail surfaces) was conducted 
in Germany in 1940 (reference 1) . The results of these tests showed 
the jet to cause (1) an increase of as much as 100 percent in the 
normal forces at small angles, (2) a shift of center of pressure to the 
rear "by an average of about 0.5 maximum body diameter, and. (3) an 
increase of drag of approximately 70 percent. Other investigations, 
both subsonic and supersonic, of jet effects upon the flow over bodies 
were conducted with the A— 4 missile at an angle of attack of 0° (refer- 
ences 2 and 3). The results of subsonic drag tests were in general 
agreement with those found in tests of the A— 5. The jet caused an 
increase of drag of as much as 80 percent. The results of the super- 
sonic tests showed a maximum decrease of drag of 18 percent. However, 
the jet appears to have been operated only at extreme over— pressure 
conditions (jet pressure as high as 60 times stream pressure) and the 
model support system causes some doubt as to the quantitative value of 
the results. 

In moat instances, test aerodynamic design of bodies housing jet 
units entails a certain degree of boat— tailing, that is, convergence, 
of the body surface as it approaches the jet exit, such that the diameters 
of the jet exit and of the exterior body surface become equal. Since 
this geometric condition would most probably favor greater jet effects 
upon the flow over the rear of the body than any other, it was chosen as 
the geometric condition to be employed in the present investigations. The 
primary purpose of the investigations was thus to determine the effects 
of an annular jet exhausting from the base of a parabolic body of revolu- 
tion upon the flow over the rear of the body. It was necessary to obtain 
first the aerodynamic characteristics of the body without the jet. 
Therefore, comprehensive force and pressure-distribution measurements 
were made of the basic jet model body. Similar, but not so exhaustive, 
auxiliary tests were conducted on a parabolic body (same body family but 
larger thickness ratio) initially employed during bench tests of small 
annular nozzles developed for use in the present investigations. All 
tests were conducted in the Langley 9— inch supersonic tunnel at a 
Mach number of 1.92. The Reynolds numbers for the tests were 2.51 X 10° 

for the jet model body and 2.47 x 10^ for the auxiliary body. 
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SYMBOLS 


A, B constants in equation of parabolic defining body shapes 

a angle of attach 

p = \]U 2 - 1 

Cp total drag coefficient / ^ )ra S \ 

\^ s max/ 


^®min 

% 

c Dp’ 

C Df 

^lam 


minimum drag coefficient 
base drag coefficient 

minimum fore drag coefficient n — 

skin— friction drag coefficient 
wave drag coefficient 

skin— friction coefficient for laminar flow on a flat 



C 


L 


c l 

^m 


^-max 

8 

€ 

l 


total lift coefficient I — \ 

y^^max ) 

weighted unit lift 

. . . /Moment about reference point 

pitching-moment coefficient 

\ ^^max^c > 


maximum body diameter 

lift density 
apex half— angle of body 
body length 
cut-off body length 
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H 

M 


M. 


des 


e 

Pa 

Pe 

Ps 


m 


Ap 

P 

*>b 

P l 


H 

Sb 

Sm 

^max 

Sw 

t 


completed body length (tip to tip) 
free— stream Mach number 

design jet Mach number based on area ratio 

angle of local surface inclination with respect to axis 
of symmetry 

radial angle 

atmospheric pressure 

static pressure of jet at jet exit 

stream pressure or pressure of ambient air 

pressure in model stilling chamber 

pressure increment 

/Ap> 


(f) 


pressure coefficient 
base pressure coefficient 
lifting pressure coefficient 


dynamic pressure 


(H 


density of' fluid 
radius of body 

Reynolds number referred to Z c 
base area 

mean cross-sectional area for body of length Z c 

maximum frontal area 

wetted area of body of length Z c 

thickness ratio 


C ^maxA 
z t ) 
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v 


volume 


V 


undisturbed stream velocity 



velocity of jet 


Cartesian coordinates 


APPARATUS AND TESTS 


Wind-tunnel and model installation .— The Langley 9— inch supersonic 
tunnel is a closed— return direct-drive type in which the pressure and 
humidity of the enclosed air may he controlled. Throughout the tests 
the quantity of water vapor in the tunnel air was kept at sufficiently 
low values so that negligible effects on the flow from condensation were 
present in the supersonic nozzle. The test Mach number is varied by 
means of interchangeable nozzle blocks forming test sections approxi- 
mately 9 inches square. A schlieren optical system provides qualitative 
visual flow observations. Eleven fine-mesh turbulence— damping screens 
are installed in the settling chamber ahead of the nozzles. 

Figure 1 shows the general installation for tests of the jet model. 
Pressure within the model stilling chamber was varied by means of 
manually controlled valves installed ahead of the juncture of the 
incoming air supply line with the flexible air supply line. Force and 
pressure— distribution measurements of the models with no jet employed 
the same model support system with the air— supply system removed. The 
scales used are self-balancing beam scales and measure three components, 
in a horizontal plane, of the total forces on the model and support 
system. 

Description of models .— All models were constructed of mild steel, 
highly polished, and excluding a special pressure-distribution model, 
were mounted on slender, hollow sting supports which, for the jet model, 
served also as an air conduit. The surface contours of the models were 
determined by revolving about its chord a parabolic arc obtained from 
the general parabolic equation 


This equation was particularly adapted for the present investigations 
since the constants A and B could be quite easily obtained for 
desired values of maximum diameter, base (or jet exit) diameter, and 
thickness ratio. (See appendix.) 


r = Ax - Bx2 


( 1 ) 
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Three separate models were constructed with a surface contour given 
by 


r = 0.1827X - 0.0l854x 2 


( 2 ) 


The designations assigned these models were : model 1-J, the basic jet 

model with two interchangeable tail sections containing jet nozzles 
of Mdes = 2.11 (nozzle 1) and M^ eg = 3-19 (nozzle 2); model 1-F, the 
model employed in the force tests; and model 1— P, a special pressure- 
distribution shell model constructed in two halves about a meridian 
plane and containing 63 pressure orifices located in one half along 
three meridians, 0°, 45°, and 90°, with 21 orifices similarly spaced 
along each meridian. 

The auxiliary model tested had a body contour given by 

r = 0.2460x - 0.02647X 2 (3) 

This model was designated model 2. 

Excluding model 1— P, the bases of all the models were hollow or 
open, as for the case of a jet exit„ Special plugs were made to fill 
the annular base openings of models 1— F and 2 flush with the body ends 
for use in tests of these models with a simulated solid or closed base. 
Photographs of models 1— J, 1— P, and 2 are shown in figure 2. Model 1— F 
has been excluded since its external appearance is no different from 
model 1— Jo 


The following table gives the pertinent geometric parameters of 
the models: 


Parameter 

Model 1 

Model 2 



7.719 

7.607 

h> in> • 

9-854 

9.293 

t 

0.09135 

0.1230 

1 , deg ........ 

IO .36 

13.83 

V 7 , cu in. . . . 

u c 

3.H3 

4.857 

Vj , cu in 

3-346 

5.080 

S v , aq in 

16.339 

20.330 

Sm, sq in. ..... . 

0.4036 

O .6385 

Sb, sq in 

O .2923 

0.3526 

S max’ 8(1 in * 

O .6365 

1.0272 

<W in * • • • 

0.9002 

1.1436 
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Development of annular nozzles .— Numerous bench tests were conducted 
to determine suitable shapes and sizes of annular nozzles that might be 
constructed in the tail section of model 1— J. Design of a theoretically 
shock— free annular supersonic nozzle contour of such small size was not 
attempted in view of the analytical complications, boundary— layer effects, 
and the difficulty of machining to the desired accuracy a curving, 
internal contour of such small radii. Nozzle 2 (M^ es = 3«19) represented 

the best attempt at construction without prohibitive surface imperfections 
of a nozzle with a curving contour to give the jet a flow direction at 
the exit similar- to that of nozzle 1 (Md es = 2.11). In spite of extreme 

care, small Imperfections in the surface contour of this nozzle could be 
detected. Because of insufficient pressure of the air— supply facility, 
conclusive bench tests of nozzle 2 could not be made. The higher ratios 
of Pe/P s (ratio of jet static pressure to pressure of ambient air) 
obtainable for nozzle 1 allowed reasonably conclusive bench tests of this 
nozzle. Diametrical surveys at the nozzle exits were conducted by means 
of a 0.010— inch total— pressure tube mounted in a micrometer— traversing 
arrangement. The total— pressure tube measured pressures on a plane 
perpendicular to the nozzle center line and just beyond the nozzle lip. 
Static pressure within the jet was measured by means of an orifice 
vented to the nozzle just inside the lip. The Mach number distribution 
across the nozzle exits was calculated from these pressures with the 
assumption of negligible effects due to the slight difference in longi- 
tudinal positions of the static— and total— pressure measurements and 
that the static pressure across the jet was constant. For the values 
of Pg/Pa of the bench tests, a conically shaped nozzle was found to 
give the most uniform distribution at the jet exit for a design Mach 
number of 2.11. Figures 3(a) and 3(^) show the results of nozzle surveys 
from bench tests and from tests conducted in a similar manner using the 
tunnel as a partly evacuated container for the model to obtain larger 
values of P e /Ps° The surveys using the tunnel as a vacuum chamber 

(p s w 0.6p a in fig. 3(a) and p a * 0.5Pa i- n fi g* 3(b)) show a marked 
improvement in the distribution for nozzle 2 and a slight lessening of 
the "hump” in the distribution curve for nozzle 1. The marked improve- 
ment in the distribution for nozzle 2 is apparently a result of the 
decrease in the pressure rise across the shock originating at the lip 
of the nozzle and reflected by the sting surface, and a decrease in the 
boundary— layer build— up caused by back pressure which in turn tends to 
eliminate compressions in the flow within the nozzle. In the bench tests 
of certain of the annular nozzles of M^ es =3 or greater (p a = p a ), 

the large pressure rise across the lip shock caused a thickening of the 
boundary layer near the lip of the outer nozzle surface and a region 
of reverse flow that extended a considerable distance away from the 
inner (sting) surface. For the cases for which reverse flow could not 
be detected, the results indicated that the large adverse pressure rise 
across the shock caused a rapid thickening of the boundary layer along 
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the sting surface ahead of the point of reflection of the shock. The 
adverse pressure gradients and the thickening boundary layers probably 
caused compressions in the flow ahead of the shock and a resulting rapid 
drop in velocity at the outer and* particularly* the inner diametrical 
stations. (See fig. 3(B).) With the decrease in external pressure 
(Ps ** 0.5Pa) the pressure rise across the lip shock decreases; therefore* 
the adverse effects would also be expected to decrease in a manner 
similar to that indicated in figure 3(B). By similar reasoning* a 
satisfactory Mach number distribution would be expected at the exit of 
nozzle 2 in the tunnel tests at M = 1.92 for which p a * 0.l4p a . 


Pressure measurements indicated that the best position for the 
orifice measuring the pressure in the model stilling chamber p^ was 

that shown in figure 1. Thermocouple measurements showed that the 
temperature of the air in the model stilling chamber varied very little 
from storage— tank air temperature. Values of the reference pressure p 

u m 

for the jet tests were measured by means of a large Bourdon type pressure 
gage. An open-tube manometer* used in conjunction with this gage* served 
as a constant check of the pressure gage and supplied values of p^ for 

pressures less than atmospheric. Figure 4 shows the calibration curves 
for each nozzle with the tunnel in operation. Although the values 
of Pq^ were intended to serve only as accurate reference pressures* 

figures 4 and 5 show that they have some quantitative value as well. 

The values of Mach number calculated from values of Pe/Pru, and presented 


in figure 4 for nozzles 1 and 2 compare favorably with the average values 
of the Mach number distributions of figures 3( a ) and 3(B)* respectively . 
In addition* figure 5 shows that the thrust of nozzle 1 obtained at two 


values of p ( 


by calculations based upon p and the Mach number 

°m 


distribution checks closely the thrust measured by strain— gage apparatus 
during the bench tests. 


Tests.— All tests were conducted through an angle— of— attack range 
of approximately ±5°° Mirrors approximately l/l6 inch square were flush- 
mounted in the bodies near the base as a part of the optical angle— of— 
attack system. Force tests and base pressure measurements of models 1— F 
and 2 were made with base open and base closed for three longitudinal 
positions of the models. These were: body base even with* 1/2 inch 

ahead of* and 1 inch ahead of the end of the sting windshield. All drag 
values were corrected for the buoyancy effect due to the difference 
between free— stream pressure and the pressure within the sting-^windshield— 
and— balance enclosing box. Radial and longitudinal surface pressure 

1 ° 

measurements were made with model 1— P at meridian intervals of 22— along 

i° 2 

every meridian from 0° to 112±- and at 180° (0° to 180° represents 
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angle— of— attack plane). With the jet in operation, the base of the model 
was 1 inch from the end of the sting windshield. The primary variables 
of the jet tests were a, p e yp 8 > and Vj ^V. For the measurements of the 

jet effects upon the pressures over the rear of the body, the tubes were 
installed as shown in the inset of figure 1. Previous investigations 
showed that the lead tubes in such an arrangement had no measurable 
effect upon the pressures over the body along a meridian l8o° opposite. 

All schlieren photographs were taken with the knife-edge horizontal. 

Precision of data .— The estimated probable errors in the aerodynamic 
quantities are included in the following table: The value of ±0.08° 

given for angle of attack is a result of error in the initial referencing 
of the model bodies with respect to stream direction. The value of ±0.01° 
Is the error that might be incurred in relative angle— of ^attack readings 
for a given test. The values for C^, Cp, and apply only to the 

results obtained from the mechanical scales. 



C D 

Cm 

M 

a 

(deg) 

R 

P 

C L 

Initial 

Relative 


± 0.0004 

±0.0004 

±0.0018 

±0.01 

±0.08 

±0.01 

±20,000 

±0.002 


Comparison of the actual ordinates of the model bodies with the values 
obtained from equations (2) and (3) showed the body dimensions to be 
accurate, with one exception, within ±0.002 inch. This exception, the 
tail section of model 1-J containing nozzle 1, had gradually increasing 
small errors in the radii of the body from a point approximately 0.3 inch 
from base rearward. The maximum error in radius (at the body base) 
amounted to ±0.008 inch. The effects of this lesser degree of boat- 
tailing will be shown in the results. The meridian planes and rotational 
angles for the radial pressure distributions were accurate within ±2°. 

The Bourdon type pressure gage for measuring Pq^ gave readings 
accurate within ±0.2 pound per square inch. 


RESULTS AND DISCUSSION 
No Jet 


Force tests .— Figure 6 shows the aerodynamic characteristics of 
model 1-F for the three longitudinal positions of the model in relation 
to the forward tip of the sting windshield. Corresponding schlieren 
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photographs for these and two additional positions are shown in figure 7* 
In like manner, the aerodynamic characteristics of model 2 and schlieren 
photographs at two longitudinal positions are shown in figures 8 and 9 , 
respectively. Excluding the zero longitudinal position, all schlieren 
photographs of figures 7 and 9 were taken at zero angle of attack. Values 
of C m in figures 6 and 8 are for moments taken about the point of 
maximum diameter. 

The results of base pressure measurements with varying longitudinal 
position and angle of attack indicated much th 9 same effects from the 
presence of the sting support and windshield and from angle of attack 
as discussed in references 4, 5, and 6 ; specifically, the base pressures 
vary appreciably with angle of attack, and the body undergoing test must 
be mounted on an extremely long, slender sting support if base pressures 
simulating free— flight values are to be obtained. The results of the 
force tests are given in table I. The base drag was given by 


The condition of the base of the bodies, open or closed, had little or 
no consistent effect upon the results except a slight increase in the 
base drag for the closed condition. The values of base drag at the zero 
longitudinal position of the models are to serve only as an indication 
of the magnitude of the fore drag and not as accurate measurements of 
the base drag since the proximity of the sting windshield to the body 
base would allow small changes of the pressure within the sting^windshield— 
and— balance enclosing box to affect the base pressures, especially for 
the closed— base condition. The effects on the lift and moment curves 
from the flow impinging upon the exposed sting at the higher angles of 
attack is shown in figures 6 and 8. These effects increase as the 
exposed sting area increases and cause marked nonlinearities in the 
curves. The longitudinal position of the models apparently had little 
effect on the minimum fore drag coefficient defined herein as 


though close examination of the schlieren photographs of figures 7 and 9 
shows a lessening of the laminar separation near the base of the body 
with increasing distance between the body base and the sting windshield. 



w 



(5) 
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Values of drag due to skin friction (laminar flow was observed over the 
entire body) were calculated for the test Reynolds number from 

S “ Cf lam(^) (6) 

These values and their approximate percentage of the fore drag are: 


Model 

C Df 

Percent of C;^ 

1— F 

0.0216 

30 

2 

.0167 

14 


References 7> 8, and 9 have pointed out independently that from 
linear theory the limiting value for the lift— curve slope of very slender 
bodies of revolution based upon is (expressed in radians) 



(7) 


and that the center— of— pressure location in relation to the nose of the 
body is 


center of pressure = 1 



(8) 


From equations (7) and (8) it follows that the slope of the pitching- 
moment curve , with moments taken about the nose of the body, is 



(9) 


The values calculated from equations (7)> (8), and (9)j expressed in 
degrees and referred to are presented subsequently and compared 

with the experimental values (in parentheses) obtained at the 1— inch 
longitudinal position. The experimental values of lift-curve and 
moment— curve slopes given in table I include support— interference 
effects and aerodynamic tares on the exposed sting. However, the 
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experimental slope values are for zero lift, and pressure measurements 
along a 1— inch length of the exposed sting from the body base have shown 
the lifting forces upon the sting to be negligible within an angle— of— 
attack range of ±2°. Furthermore, the effects of the presence of the 
sting and windshield upon the body lifting forces would be expected to 
be least at the 1— inch position. 


Model 

dC L 

da 

c.p. 

(diam from nose) 

dCm 

about nose 

da 


0.0160 

3.26 

-0.00610 

1-F 

(.0290) 

(1.28) 

(-.00431) 


.0120 

5.49 

-.00975 

2 

(.0338) 

(1.39) 

(-.00705) 


All of the theoretical values are relatively poor predictions of the 
experimental results. Part of the failure of equations (7)* (8), and (9) 
to predict values in reasonable agreement with experimental values is 
probably caused by the use of the geometrical value of base area. 
Excluding the very slenderest of bodies and the case where the base area 
approached the frontal area, the flow over the rear of parabolic bodies 
of revolution or similar body shapes can hardly be expected to take 
place without some form of separation. Eeference 4 has shown that the 
calculated pressures over the rear of a body of revolution with boat- 
tailing, as given by the method of characteristics, are in excellent 
agreement with experimental pressures if the pressure calculations are 
performed along the streamline of separated boundary layer. This would 
seem to indicate that the geometrical value of base area in equations (7)* 
(8), and (9) should be replaced by an area determined by the diameter 
between the separated streamlines at the body base. Measurements of 
this "diameter of separation" were made from enlarged schlieren photo- 
graphs of models 1— F and 2, each at the 1— inch longitudinal position 
and with open base. An area of 0.347 square inch was obtained for 
model 1-F and an area of 0.47^ square inch for model 2. The values 
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calculated from equations (7 )j (8) , and (9) and referred to these areas 
are presented in the following table and compared with the experimental 
values (in parentheses): 


Model 

ac L 

da 

CoP. 

(diam from nose) 

dQm 

about nose 

da 

1— F 

0.0190 

1.40 

-0.00310 


(.0290) 

( 1 . 28 ) 

(-.00431) 

2 

.0161 

2.31 

-.00559 


(• 0338 ) 

( 1 . 39 ) 

(-.00705) 


These values are an improvement upon the previous theoretical values, 
but the lift predictions are still rather poor, which in turn affect the 
pitching-moment values. 

Pressure distributions .— The results of radial pressure-distribution 
measurements are presented In figure 10 for model 1— P and in figure 11 
for model 2. Longitudinal pressure distributions are presented in 
figures 12 and 13 for model 1— P and in figure 14 for model 2. Although 
the results for model 2 are secondary to those for model 1— P, they tend 
to indicate that certain phenomena observed in the pressure distributions 
of both bodies apparently hold for slender pointed bodies of revolution 
in general. Some of these phenomena have been previously observed in 
reference 6 and, in particular, in the investigation of a typical super- 
sonic aircraft fuselage of reference 10. First, figures 10, 11, 13(b), 
14(c), and 15 show that the pressures along the 90° meridian at a = 0° 
do not remain relatively unchanged with angle of attack - , a simplifying 
assumption generally employed in existing approximate theories for 
computing the aerodynamic characteristics of conical bodies and pointed 
bodies of revolution. In fact, at certain horizontal stations the 
pressure at the 90° meridian varies as much or more than at any other 
meridian. Second, the radial pressure distribution at any longitudinal 
station varied appreciably from the usually assumed cosine distribution, 
especially so for the low-pressure half of the body at longitudinal 
stations ahead of the maximum thickness and for the entire circumference 
at stations behind the maximum thickness. As noted in reference 10, 
the radial pressures at stations after the maximum thickness behave in 
much the same manner as observed in tests of slender cylinders in yaw 
(reference 11). Third, for longitudinal stations ahead of the maximum 
thickness there appears to be a radial point for each longitudinal 
station at which the pressure remains almost constant with angle of 
attack - . (See fig. 15, in particular.) The locus of these radial points 
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does not follow a meridian but tends to move away from the 90° meridian 
in the direction of the high-pressure half of the body as the distance 
from the nose increases. For model 1— P, this shift was from 9 ~ 80° 
at station 0.088 to 9 ~ 48° at station 0.606. For model 2, the shift 
was from 0 £ 78° at station 0.283 to 0 & 6l° at station 0.573* The 
fourth phenomenon observed was the consistent appearance of the "hump," 
not predicted by theory, in the longitudinal pressure— distribution 
curves. For model 1— P, this characteristic occurred near the 0.3 station 
and for the a = 0° condition amounted to a noticeable discontinuity 
in the curve. For model 2, it occurred in the vicinity of the 0.5 
to 0.6 stations. This phenomenon has been present in the results of 
other tests of slender pointed bodies of revolution (references 6, 10, 
and 12, for example) and, excluding cones, is apparently characteristic 
of slender pointed bodies of revolution in general. All of the above 
phenomena are undoubtedly associated with the complicated nature of the 
flow and the viscous effects that exist in the flow over inclined 
slender pointed bodies of revolution. Though cross velocities in the 
vicinity of the 90° meridian would not be expected to affect the lift, 
their inclusion would, nevertheless, be expected to reduce the pressures 
at the 90° meridian, possibly of the magnitude observed in the experi- 
mental results. Also, the fact that when the experimental longitudinal 
pressure gradients in the vicinity of the 90° meridian are found to be 
relatively large, the experimental tangential pressure gradients are 
found to be of the same order of magnitude. A tangential gradient of 
such magnitude would be expected to have important bearing upon 
separation effects. 

Figure 16 illustrates the method by which the pressure coefficient 
at any point on the body is converted to the lifting pressure 
coefficient P^. The equation, including second— order terms, is 


All values of P for model 1— P were converted to P^ by means of 
equation (10). With as the reference area, the total lift 

coefficient would be given by (see top of fig. 17) 


P2 = P(cos 0 cos q cos a + sin q sin a) 


( 10 ) 



cosq 


(ID 
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Plots of P^cos tj against 6 were made for each horizontal station 

and graphically integrated from 0 to 2jt. This gave what might be 
termed the lift density 6 at each station. Therefore, the weighted 
unit lift for each station x is 


( c i) x = Vx ( 12 ) 


From the linear first approximation theory of reference 8 and the first 
approximation theory of reference 13, a solution for the lift distribu- 
tion over the body from equation (7) would apparently give an acceptable 
first— order prediction. In the dimensions of equation (12) and for 
small values of a, this solution may be expressed as 

( c j) x - 1,OT < r >*(£) <«> 

x 


In figure 17 the values of c^ from equations (12) and (13) are plotted 
against horizontal station, in inches, for values of a of 2.50° 
and 5. 00°. This gives a graphical representation of the lift distri- 
bution over the body. Values of Cp and C m determined by integration 
of the experimental curves and the theoretical values in parentheses 
are presented in the following table: 


a 

°L 

C 

m 

c.p. 

(&eg) 

about nose 

(diam from nose) 

2.50 

O0O7I7 

-0.00733 

0.88 

(. 0400 ) 

(-.01525) 

(3.26) 

5.00 

.1472 

-.0230 

1.34 

(.0800) 

(.03050) 

(3.26) 


Since the tests of model 1— P most nearly duplicate in tunnel 
position the tests of model 1-F at the 1— inch longitudinal sting station, 
a reasonable check of the above pressure— distribution results should be 
realized in a comparison with the force data results for the 1— inch sting 
station (open base). In an effort to eliminate as much as possible the 
effects of the sting support, values based upon the slopes of the 
curves at zero lift were determined from table I. These values are 
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included in the following table. Of course., the use of constant slopes 
will mask the effects of separation at the rear of the body. 


a 

(deg) 

°L 

Cm 

about nose 

c.p. 

(diams from nose) 

Source 

2.50 

O0O725 

-0.0108 

1.28 

Constant Slope 

5.00 

.1450 

-.0216 

1.28 

Constant Slope 

2.50 

.0700 

-.00766 

.94 

Average— curve value 

5.00 

.1450 

-.0146 

.86 

Average— curve value 


In an effort to obtain better agreement between the results of force 
tests and pressure measurements at a = 5*00° (at the 1— inch longi- 
tudinal sting station) model 1 -F was tested as shown in figure 18. 
The results are given in figure 19 and tabulated in the following 
table : 


a 

°L 

Cm 

c.p. 

(deg) 

about nose 

(dtam from nose) 

2.50 

0.0860 

-0.0246 

2.45 

5.00 

.1831 

-.0532 

2.49 


These values show no close agreement with the values obtained from 
model 1 — P, but it is interesting to note that they agree closely with 
the results of the tests of model 1-F, open base., at the zero longi- 
tudinal position, thus indicating an appreciable effect from the flush 
sting windshield arrangement. As might be expected, the only close 
check of the values obtained through integration of the pressure 
distributions lies in a comparison at a = 2.50° with the average- 
curve values from tests of model 1-F at the 1— inch longitudinal sting 
station. On this basis, the results of the pressure-distribution 
integrations may be considered reliable. The shift of center of 
pressure with angle of attack, as determined from the pressure distri- 
butions, would correspond to effects that might be expected from 
separation of the flow at the rear of the body. 

The results of the pressure— distribution investigations give some 
insight into the causes of the higher experimental values of dC^/da, 
as compared with theoretical values, noted in these and other tests of 
slender pointed bodies of revolution. Figure 17 shows that, for body 
stations behind the station at which the theoretical local lift would 
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be a maximum, the experimental values of local lift exhibit a marked, 
increase over the theoretical values. The station of initial gain in 
the experimental local lift also appears to be in the region where the 
experimental radial pressures begin to deviate appreciably from a cosine 
distribution. It is seen that the observed non— cosine distributions of 
radial pressures act in such a way as to give more lift over the body 
than would cosine distributions. Figure 17 also shows the expected 
reduction of anti— lift forces by separation in the region over the rear 
of the body where recompression would otherwise take place. Reference 14 
has shown that the effect of the boundary layer alone (no separation) 
is to increase the ordinary linear theory result by an amount dependent 
upon the displacement thickness. The sum of these effeqts could 
appreciably increase the nonviscous theoretical value for a lift— curve 
slope of 2 (per radian) . 

Analysis of theoretical methods for prediction of longitudinal 
pressure distribution .— The equations for several methods for predicting 
the pressure distribution over slender pointed bodies of revolution were 
calculated in terms of equation (1). Reference 9 gives a solution 
termed the "rigorous linearized first-order solution" for an extremely 
sharp— nose body of revolution for M > \/2L This yields 


P « — 2(a 2 — 6ABx + 6B 2 x 2 )log(A — Bx) 


(14) 


which was obtained in reference 9 from a previously derived equation for 
the pressure coefficient for compressible flow, given herein as 


P = -A 2 + l6ABx - 22B 2 x 2 + 2(A2 - 6ABx + 6B 2 x2) 



2 - logOA - pBx)] 


-I- 


0(a3 - 3A2bx + 3AB 2 x2 - b3x3) 


(15) 


By the method of reference 15, but with axes at the nose of the body, 
the general equation was found to be 


P = 


2 < (6ab 



P 2 r 2 + £a 2 — 6ABx + 3®^(2x 2 + 



(16) 
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Reference 16 gives a solution for P that is identical with equation (16) 
except for one additional term. This solution is 

P = (value from Equation (16)) — (A^ — 4ABx + 4 B^x2) ( 17) 

Reference 12 has given a solution "based upon the small— di sturhance theory 
and requiring a step— by— step numerical integration. For simplicity this 
method will he expressed herein merely as 


P = 


]T|hr,Mj] 

i=l 


(18) 


where i and n represent the initial and nth integration station, 
respectively. A 25— point calculation was made. The final method 
employed was the characteristics method of references 17 and 18 with the 
assumption of potential flow. 

Equations 14 to 18 were applied to the body shapes of models 1 
and 2. The characteristics method was applied to model 1 only. The 
results of these pressure-distribution calculations are presented in 
figure 20. The corresponding wave drag coefficients Cp^ from inte- 
gration of the pressure curves are given in the following table: 


Method 

Model 1 

Model 2 

Equation (14) 
Equation (15) 
Equation (16) 
Equation (17) 
Equation (l8) 
Characteristics 
method 

0.0253 

.0784 

.0822 

.0740 

.0746 

.0624 

0.0356 

.1498 

.1548 

.1302 

.1309 


Integration of a curve determined by the experimental points of 
figure 20 gives for model 1 a pressure drag coefficient of 0.0542 and 
for model 2 (very approximate due to lack of sufficient points) a drag 
of 0.11. It should be noted that no correction has been applied to the 
experimental points since surveys have shown negligible variation in 
the static pressure distribution in the wind-tunnel test section. 
Therefore, any correction applied to the pressure drag would be 
negligible. The effects of separation upon the experimental pressure 
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drag coefficient would not oppose the indication that all of these 
theories predict too large a wave drag. As might he expected., the 
characteristics method shows test agreement with experiment. Though it 
gives only a fair prediction* equation (17)* the Lighthill solution* 
is the best of the approximate theories and gave a slightly better 
prediction of the experimental results than did the laborious small- 
disturbance method of equation (18). Indications are that equation ( Ik ) 
should not be applied. 

If the values of the experimental pressure drag are assumed to 
approach the wave drag* the addition of the laminar shin— friction drag 
should give a value that checks closely the measured fore drag. The 
comparison is given in the following table. The corresponding values 
of the fore drags are from table I* 1— inch position* open base. 


Model 

Experimental pressure 

Fore drag from 

drag plus Cj ^ 

force tests 

1 

0.0758 

0.0743 

2 

.1267 

.1104 


With Jet 

Figure 21 presents schlieren photographs of the jet model with 
and without tubes to the pressure orifices installed as shown in the 
upper left— half corner of figure 1. As previously stated* pressure 
measurements with no jet throughout the angle-of-^attack range showed 
no effect upon the body pressures from the presence of or disturbances 
produced by the pressure lead tubes on the side of the body opposite 
the pressure orifices. The surveys and calibrations of the two jet 
nozzles indicated reasonable values of the Mach numbers for the two 
jet nozzles to be approximately 2.10 and 3°05<> For a free— stream Mach 
number of 1.92* these values would represent ratios of jet velocity to 
stream velocity Vj^V of approximately 1.05 and 1.24* respectively. 

Figure 22 shows the pressure change at each orifice location due 
to jet effects with varying jet pressure and angle of attack. Also 
included are the hysteresis effects (at the a = 0° and a = 2.50°, 
upper surface, for both velocity ratios) dependent upon whether the 
particular test was made with increasing or decreasing jet pressures. 
For both velocity ratios the major effects upon the pressures over the 
rear of the body occurred at the a = 0° and a = 2.50° upper-surface 
conditions and were confined to the orifices nearest the jet exit. At 
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these a conditions and a velocity ratio of 1.05* the body pressures 
increased positively as the jet pressure approached and exceeded stream 
pressure , the greatest change occurring immediately after the over— 


and a velocity ratio of 1.24, the body pressures showed a very slight 


at the over— pressure conditions. Figure 23 gives the percentage change 
in body pressures due to the jet at a = 0° and a = 2.50°, upper 
surface , for several values of P e /Ps* At the top of the figure the 
differences in the basic pressure distributions over the* rear of model 1— P 
and the two nozzle tail sections (no jet) of model 1— J are presented for 
a = 0°. These differences in pressures appear to be the effects of 
previously mentioned very small measured differences in body ordinates. 

The zero— percent datum lines of the plots of jet effects represent the 
pressures with no jet. The major effects of the jet upon the body 
pressures are confined to approximately 5 percent of the body length 
(from the base) for a velocity ratio of 1.24 and to approximately 10 per- 
cent of the body length for a velocity ratio of 1.05. For similar 
pressure ratios p e /p Q the effect of the jet upon the body pressures in 
these regions is much greater for the lower velocity ratio. At a = 2.50°* 
upper surface, there is a positive increase in jet effects over the 
a = 0° condition.. This might be expected in view of the greater separa- 
tion of the flow from the upper surface at angle of attach- that would 
favor pressure effects from the jet to be felt farther forward along the 
body and to a greater degree. 

Typical schlieren photographs of the jet in operation at a = 0° 
are shown in figure 24 for the two velocity ratios and, whenever possible, 
for similar pressure ratios. Photographs at the lower velocity ratio 
were taken with the lead tubes installed; therefore, for comparison with 
the photographs at the higher velocity ratio which were taken without the 
tubes installed, the simple pattern of disturbances present in the no- 
jet photograph should be ignored. 

As the supersonic jet begins to flow, there is a noticeable decrease 
in thickness of the laminar separated region at the rear of the body. 

This is particularly true at the higher velocity ratio. Up to a value 
of Pe/Pg equal to approximately 0.5 the shock pattern within the jet 
and at the nozzle lips is much the same for both velocity ratios. A very 
prominent lambda shock is noted at the jet outer boundary immediately 
rearward of the nozzle lip. The shock pattern within the jet follows 
closely the expected phenomena for under— pressure jets calling for the 
presence of a shock to offset the higher pressure outside the jet 

boundary. From = 0.5 to 1.0 the lambda shock at the jet outer 


pressure condition ^ > 1 was reached 


At the same a conditions 


decrease at the under— pressure conditions 



1 and a slight increase 
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■boundary tends toward a plain shock whose apparent point of origin at the 
Jet boundary lies slightly downstream of the nozzle lip. The shock 
pattern within the Jet continues along the pattern for under— pressure 
Jets except that two shocks are now observed within the Jet of higher 
velocity ratio. It is possible that this may be due to slight imperfec- 
tions in the nozzle surface; but, in view of the fact that a similar 
phenomenon* though not as strong* may be observed at slightly 
higher pressure ratios at the lower velocity ratio* it appears more 
probable that this phenomenon is characteristic of the flow from 

annular nozzles of this type. At a value of ^2- w 1 the Jet boundaries 

f s 

are parallel to the axis of symmetry of the Jet. Parallel flow would 

not be expected to occur exactly at = 1 because the stream pressure 

would be somewhat different from the pressure behind the shock emanating 
from the vicinity of the nozzle lip. Also* the presence of the shocks 

within the parallel Jet at ^ ~ 1 may be attributed to the necessary 

Dg 

change in flow inclination from a direction that is away from the axis 
of symmetry immediately ahead of the exit at the nozzle lip to a 
direction parallel to the axis of symmetry beyond the exit. At values 


of > 1 the separation over the rear of the body increases* this 
-^s 

effect being more pronounced at the lower velocity ratio. An expansion 
of the jet flow occurs at the nozzle lip for the over-pressure 

conditions ^- > 1 and is followed by shocks within the jet. 
rs 


(See 


fig. 25.) Almost identical phenomena have been observed in reference 19 
in the flow around the trailing edge of a flapped airfoil from the high- 
pressure side. The flow apparently overexpands and the shock is 
necessary to raise the pressure to satisfy the jet— boundary conditions. 
With increasing over— pressures of the jet, the shock from the outer 
boundary of the jet near the jet exit reverts more and more to a lambda 
shock at the higher velocity ratio, whereas the corresponding shock at 
the lower velocity ratio continues as a clearly defined plain shock. 
Indications from the pressure measurements and the schlieren photographs 
are that the induced velocities imparted to the dead air in the region 
where the flow is separated from the rear of the body are greater for 
the case of the higher velocity ratio. Apparently, these induction 
effects counteract the back— pressure effects over the body from the 
related internal (jet) and external flow characteristics. 
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The maximum effects of the jet upon the fore drag of the body 
(a = 0°) were determined from the body pressures with the jet in opera- 
tion. The results are tabulated in the following table as percentage 
change in Cp^, : 


The maximum effect of the jet upon the lift and pitching moment 
occurred at a = 2 . 50 ° and, as in the case for the fore drag, at the 
lower velocity ratio, 1.05 j and highest pressure ratio, 3*5« An 
approximation of this effect is entered on the curve representing lift 
distribution in figure 17. The results of such a change in lift 
distribution would be a 13 . 7 - percent decrease in the over-all lift of 
the body and a destabilizing shift of center of pressure from 0.88 diameter 
rearward of the nose to a point 0.74 diameter ahead of the nose. The 
foregoing determination of jet effects has assumed the interference 
effects of the sting and windshield on the pressure measurements to be 
the same for the jet— on and jet-off conditions. However, the differences 
in the effects are expected to be small. Thus, the foregoing values 
should give an insight into the magnitude of some of the effects that 
might be expected from an exhausting annular jet. 

While the present tests to determine jet effects for a given design 
condition cover the range of pressure ratios P e /Pg that might be 
expected for actual flight conditions, they are limited in the range of 
velocity ratios V j/V. Future investigations of this sort should attempt 
to achieve velocity ratios of approximately 2 in an effort to duplicate 
flight conditions. There is a definite need for fundamental investiga- 
tions of supersonic mixing zones such as stem from jet exits. 


Supersonic investigations at a Mach number of 1.92 of a parabolic 
body of revolution with and without an annular jet exhausting from the 
base indicate the following conclusions: 



0.4 0.34 

1.0 -.21 0.72 

1.5 - 1.01 0 

1.8 - 1.62 -.43 

3-5 -3.31 


CONCLUSIONS 
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For the case with jet inoperative, 

1 . The condition of the t»ase of the tody, hollow or closed, had 
little and no consistent effect upon the aerodynamic characteristics of 
the tody. 

2 . The simplified linearized solutions for lift-curve slope, 
pitching-moment-curve slope, and center— of— pressure location gave poor 
predictions of the experimental results. An improvement upon the 
theoretical values was obtained ty substituting an n area of separation" 
at the tody base for the base area in the equations. 

3 . Experimental radial pressure distributions showed a marked 
deviation of lifting pressures from the theoretical cosine distribution, 
and, contrary to the simplifying assumptions of approximate theories, 
the pressures on the sides of the body (90° from angle-of-attack plane) 
varied appreciably with angle of attack. It is quite possible that 
these discrepancies are a result of separation effects in the thick 
laminar boundary layer and of the failure of the theories to include 
effects of cross velocities which may be important. 

4 . The method of characteristics for axial symmetry gave a 
reasonable over-all prediction of the actual pressure distribution over 
the body. However, a "hump," not predicted by theory, was found in the 
experimental longitudinal pressure— distribution curve at forward body 
stations. This phenomenon appears to be characteristic of slender 
bodies of revolution in general. 

5 . The Lighthill solution appears to be the best of the linearized 
solutions investigated for prediction of pressure distribution over 
slender bodies of revolution. 

For the case with Jet in operation, 

6. Pressures over the rear of the body showed little effect from 
the jet until the jet pressure approached and exceeded the value for 
parallel flow of the jet. 

7 * The effects from the jet upon the body pressures were reduced 
as the ratio of jet velocity to free— stream velocity increased. 

8. Maximum effects of the jet were obtained at the lower ratio 
of jet velocity to stream velocity, 1.05* and the highest ratio of 
jet pressure to stream pressure, 3 * 5 * These effects amounted to a 
3.3— percent reduction in fore drag at 0° angle of attack and, 
at 2. 50° angle of attack, a 13. 7— percent reduction in lift and a 
destabilizing shift of center of pressure from 0.88 body diameter 
rearward of the nose to approximately 0.74 diameter ahead of the nose. 


CONFIDENTIAL 


24 


CONFIDENTIAL 


NACA EM L9K09 


9. The greatest effects of the jet upon the lift and pitching 
moment occurred at 2.50° angle of attack- and almost completely dis- 
appeared as the angle of attack - was increased to 5»00°. 

10. Indications are that, for higher ratios of jet velocity to stream 
velocity than achieved in the present investigations, the induction 
effects of the jet upon the flow over the rear of the tody would become 
more important than tack— pressure effects. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX 


The general equation for the shape of the todies is given as 


r = Ax — Bx^ 


(Al) 


The constants A and B may he determined simply, as follows: 
At maximum thickness 


— - = A — 2Bx = 0 
dx 


therefore. 



or, in nondimensional form 


r max _ A^ 

l 4b i 


(A2) 


Similarly, 


r tase(or jet exit) 


Al - B l 2 


whence 


A — laae 
B = i — 

l 


(A3) 


r max-> r base-» and- thickness ratio are assigned fixed values, 

the constants A and B are readily obtained by combining equations 
(A2) and (A3). Examination of the basic equation (Al) shows that the 
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constant A is dimensionless and is equal to twice the thickness ratio. 
However, the constant B is not dimensionless since it hears the 
following relation: 


B 


A2 


hr 


max 


which has the dimension 


1 

unit of length 


Therefore * any calculations employing equation (Al) with the numerical 
values replacing the constants A and B, such as computations of 
pressure distributions or drag., must he carried through with the same 
dimension units (ft or in,) used in calculating the values of A and B. 

The equations given below apply to the family of bodies determined 
from equation (Al) 


Volume 



ABZ 2 * + B g Z 5 

2 5 


(Ak) 


Mean area 


A 2 z 2 _ ABz3 + b2z^ 
3 2 + 5 


(A5) 


Wetted surface area 


n |az 



(a6) 
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TABLE I.- SUMMARY OF FORCE TESTS RESULTS FOR 


MODELS 1-F AND 2 WITHOUT JET 


Model 

Distance from 
"body base to 
forward end of 
sting windshield 
(in.) 

Base 

condition 

{-) 

\da y L=0 

at ° ut Vi 

CD min 

No>0° 

c.p, 

(diam from nose) 


0 

Open 

0.0334 

0.0121 

0.0899 

0.0142 

2.37 


0 

Closed 

.0330 

.0124 

.0879 

.0156 

2.25 

th 

0.5 

Open 

.0324 

.0126 

.0940 

.0214 

2.14 

1-F 

0.5 

Closed 

.0323 

.0102 

.0955 

.0233 

2.77 


1.0 

Open 

.0290 

.0142 

.0960 

.0217 

1.28 


1.0 

Closed 

.0292 

.0140 

.0960 

• 0235 

1.45 


0 

Open 

.0392 

.0112 

.1370 

.0217 

2.16 


0 

Closed 

.0394 

.0110 

• 1373 

.0220 

2.21 

p 

0.5 

Open 

.0427 

.0105 

.1400 

.0241 

2.43 

C. 

0.5 

Closed 

.0429 

.0105 

.1400 

.0252 

2.43 


1.0 

Open 

.0338 

.0136 

.1410 

.0306 

1.39 


1.0 

Closed 

.0343 

.0130 

.1407 

.0309 

1.54 
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Figure 1.- Jet model tunnel installation. 
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(a) Models 2 and 1— J. 



(b) Side view of model 1— P. 



(c) 


Auxiliary view of tube exits, model 1— P. 


Figure 2.- Models. 
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(a) Nozzle 1, M des = 2.11. 
Figure 3 -- Surveys of nozzle exits. 
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(b) Nozzle 2, M des = 3*19* 
Figure 3.- Concluded. 
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Figure lj-.- Calibration of jet nozzles in tunnel. 


T/irustj pounds. 


3 6 


NACA EM L9K09 


CONFIDENTIAL 



Figure 5*- Thrust of nozzle 1. 
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(a) Base to windshield is 0 inches. 

Figure 6.- Aerodynamic characteristics of model 1-F. 
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(b) Base to windshield is 0.5 inches. 
Figure 6.- Continued. 
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(c) Base to windshield is 1.0 inches. 
Figure 6.- Concluded. 
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(a) 0 inch. 



(b) 0.5 inch. 



(d) 1.5 inches. 



(c) 1 inch. 



(e) 2 inches. 



L-630li2 


Figure 7* - Schlieren photographs shoving the effect upon the flow over 
the rear of model 1-F with varying distance between body base and tip 
of sting windshield. 
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(a) Base to windshield is 0 inches. 

Figure 8.- Aerodynamic characteristics of model 2. 
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(b) Base to windshield is 0.5 inches. 
Figure 8.- Continued. 


NACA EM L9K09 


^5 



(c) Base to windshield is 1.0 inches. 


Figure 8.- Concluded. 
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(a) 0 inch. 



(b) 0.5 inch. 


Figure 9-- Schlieren photographs shoving the effect upon the flow over 
the rear of model 2 with varying distance between body base and tip of 
sting windshield. 
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Figure 10.- Radial pressure distributions for model 1-P. 
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Figure 10.- Concluded. 
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(a) Body stations ahead of maximum thickness station. 
Figure 11.- Radial pressure distributions for model 2. 
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(b) Body stations behind maximum thickness station. 


Figure 11.- Concluded. 
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&= 90 ° 





Figure 12.- Variation of longitudinal pressure distribution with angle 
of attack at two meridians of model 1-P. 
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Figure 13.- Longitudinal pressure distribution of model 1-P at a = 
and variation of 90° meridian pressures with angle of attack. 
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Figure 14.- Longitudinal pressure distributions for model 2. 
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8 = 270 ° 
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Figure 15 .- Radial pressure distributions for model 2. (Examples of 
radial point of zero pressure variation with angle of attack. ) 
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Figure 1 6 .- Conversion of normal pressures to lifting pressures. 
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Figure 17- - Graphical representation of the lift distribution over 

model 1-P. 
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Figure 18.- Model 1-F with flush sting windshield. 
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Figure 19 •- Aerodynamic characteristics of model 1-F, surface of sting 

windshield flush with body base. 
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Figure 20.- Comparison of theoretical and experimental pressure distribu- 
tions for a = 0°. (Numbers in squares indicate equation numbers in 
text. ) 


NACA PM L9K09 


63 


Without pressure 
lead, tubes 


CONFIDENTIAL 


Pressure lead 
tubes Installed 



Pe 

(c) = 1.71. 

-^8 


^NACA. 

l-63014j- 


(c) — = 1.70. 

P s 


Figure 21.- Schlieren photographs of phenomena at jet exit with and vith- 
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out pressure lead tubes installed; = 1.2h. 
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Figure 22.- Effects of jet upon pressures at rear body stations with varying 

angle of attack and jet pressures. 
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Figure 22.- Concluded. 
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Figure 23-- Effects of jet upon rear body pressures. 
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(a) No jet. 



(a) No jet. 



Figure 24.- Schlieren photographs of jet in operation at two ratios of 
jet velocity to free- stream velocity and with varying ratios of jet 
pressure to free-stream pressure. 
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Figure 2k. - Continued. 
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Figure 2b. - Concluded. 
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(a) Schlieren photograph of flow phenomena. 



(h) Sketch of flow phenomena. 
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Figure 25.— Phenomena in region of jet exit* = 1.05* — = 3*1°* 
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